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In underdoped cuprates, is thought to be determined by Cooper pair phase 
fluctuations because of the small superfluid density Us- Experimentally, Tc 
is found to scale with Ug. The fluctuation- exchange approximation (FLEX) 
in its standard form fails to predict this behavior of Tc since it does not 
include phase fluctuations. We therefore extend the FLEX to include them 
self consistently. We present results for Tc[ns,x], where x is the doping. 
PACS numbers: 74.72-h, 74.20.Mn, 74.25. Jb 



One of the striking properties of underdoped cuprates is the so-called 
strong pseudogap above Tc.tl It appears below a characteristic temperature 
T* in the underdoped regime and T* approaches the superconducting transi- 
tion temperature Tc near optimal doping, where Tc is maximal. The pseudo- 
gap evolves smoothly into the superconducting gap and has the same d-wave 
symmetry. A further characteristic property of underdoped cuprates is the 
low superfluid density Ug. Furthermore, T^ oc ns{T = {))/m* is observed 
experimentallyi and there are indications of finite phase stiffness being pro- 
portional to Us on short length and time scales in the pseudogap regime 
above Tc.i One interpretation is that at T* Cooper pair formation occurs 
and that at Tc these Cooper pairs become phase coherent. In accordance 
with this physical picture we extend FLEX to include phase fluctuations. 



Theory 



As a microscopic model for high-Tc superconductors we use the single- 
band Hubbard model in 2D, with the nearest-neighbor hopping element t and 
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the local Coulomb repulsion U . The fluctuation exchange approximationQ 
(FLEX) is used to treat the superconducting state in the standard Nambu- 
Eliashberg formalism with the spin fluctuations of the FLEX as the pairing 
interaction.!' i From the self-consistent solution of the Dyson equation one 
can determine T* as the highest temperature for which the off-diagonal self 
energy is non- vanishing. Usually this is the criterion for the critical tempera- 
ture Tc, but this is only true if one can neglect the role of phase fluctuations. 
However, when the superfluid density is small, phase fluctuations will destroy 
phase coherence, while the charge carriers still form local Cooper pairs, which 
is signaled by a non zero amplitude of the superconducting order parameter. 

To extend the FLEX to include phase fluctuations self-consistently, one 
has to know their self-energy contribution. In order to derive it we start with 
a phenomenological model of tight-binding electrons and an effective nearest- 
neighbor pairing interaction for electrons with opposite spins, leading to a 
superconducting c?-wave order parameter in the mean- field approximation. 
We begin by writing down the action S for the tight-binding electrons where 
the interaction term is decoupled by a Hubbard-Stratonovich transformation, 
introducing a decoupling field Ajj , which turns out to be the superconducting 
order parameter: 

S[^\ CD; A*, A] = f drlY. ^-(^- " ^)^- + * E 
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where represents the electron field at lattice site i with spin a. V is 
the pairing interaction and t is the hopping matrix element between nearest- 
neighbor sites. The doping is controlled by the chemical potential /i. We 
now assume the local superconducting order parameter Aj^ (r) to have a time 
and translational invariant amplitude A^j with d-wave symmetry but a fluc- 
tuating phase suppressing superconducting order in the pseudogap regime 
below T* . We perform the following transformation to decouple the phase 
and amplitude degrees of freedom of the order parameter J3 

= e-5^»W , A,, = AO. e-^I^^M+'^.Wl . 

This leads to an expression for the action where the phase is directly coupled 
to the fermions by the elementary vertices shown in figure ||. After integrat- 
ing out the fermionic degree of freedom, a loop expansion up to second order 
in the phase fields ipiir) leads to the effective action for the phases, taking 
the form of a time-dependent XY model on a cubic lattice. From the effective 
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Fig. 1. In (a) the elementary vertices ■ and • of our theory are shown. 
The solid line represent an in or outgoing electron and the wiggly the phase 
field. In (h) we show the diagrams contributing in lowest order to the effec- 
tive action 5e//. for the phases, here the solid line is the Namhu Green 
function. 

action we determine the propagator of the phase fluctuations: 

W"^'{q^iu^) = — . (1) 

-K{q,iun) {ii^ny + '^T" 2(2 - cosg^ - cosg^) 

where K{q,iun) is the density-density correlation function and the dynam- 
ical phase stiffness is given by ns{q,ivn) /m* = T + {q,oj) / e^ . Here T 
is the expectation value of the inverse band-mass operator and I!"'-' is the 
current-current correlation function calculated in the mean-field approxima- 
tion. The superfluid density is the density of phase coherent Cooper pairs 
only, which are responsible for the Meissner effect. If no phase fluctuations 
are present it will take the value of the local Cooper-pair density Ucp. Since 
in the London theory Ug is related to the London penetration depth by 
Xi = Y^m*c?/(47r€?n^ we can determine Tc as the temperature where the 
Meissner effect and therefore the superfluid density ns(g = 0, iz/„ = 0) refer- 
ing to phase coherent Cooper pairs vanishes. This can be seen as the loss 
of rigidity against long-range phase fluctuations and is consistent with the 
transition from quasi long-range to short-range order in 2D.i For T > Tc 
the dynamical phase stiffness will be non-zero for certain q > qo and v > i'q 
and will result in superconducting order on finite length and time scales and 
will finally vanish for all momenta and frequencies at T*, where Ucp goes 
to zero.B From the generating functional for the fermions we determine the 
self-energy corrections due to phase fluctuations in lowest order. Combining 
this self-energy contribution with the contribution of standard FLEX and 
taking the renormalized propagator of the phase fluctuations, which means 
that one has to calculate ns{q,iyn) and K{q,iun) with renormalized Green 
functions, gives a set of self-consistency equation for the renormalized Green 
functions. In this way we get an extension of the usual FLEX equations with 
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Fig. 2. Dyson equation for the extended FLEX approximation. The straight 
double lines are renormalized Nambu Green function, whereas the single lines 
are bare ones. The zig-zag line in the second graph is the effective interaction 
of the FLEX, and the wiggly lines denote the propagator H'^'^ of the phase 
fluctuations. 

the Dyson equation shown in figure |^. 

Results 

In figure ^ we show results for the critical temperature T* at which 
Cooper pairs are formed and for Tc where the phase of the Cooper pairs 
become coherent. T* is determined within standard FLEX as the highest 
temperature for which the off-diagonal self energy is non-zero. Tc is calcu- 
lated using the Kosterlitz-Thouless criterion0 7r/2 h^dns{Tc)/^m* = ksTc 
for the temperature where phase coherence vanishes. Here ns(T)/m* is cal- 
culated within FLEX approximation and d is half of the c-axis lattice con- 
stant. Note the Kosterlitz-Thouless criterion neglects the coupling between 
amplitude and phase fluctuations, which is expected to be important for the 
overdoped regime. Also we show FLEX results for n^/m* extrapolated to 
T = 0, indicating Tc oc nJT = 0)/m* in the underdoped regime. We flnd 
Tc = 1.691/i^dns/4m* /cbJ^j In accordance with experimental findings we ex- 
pect that our results for ns{T = 0)/m* and Tc should decrease more rapidly 
for doping a; — > due to the vicinity of the antiferromagnetic transition 
which is not properly described in the FLEX. It is interesting to note, that 
Tc T* in the overdoped regime as expected for increasing Cooper pair 
density, although using the Kosterlitz-Thouless criterion neglects coupling 
of amplitude and phase fluctuations. We obtain Tc ~ T*, since ng/m* in- 
creases rapidly below T* as compared to the underdoped regime. The result 
for ns{T = 0)/m* in the overdoped regime suggests that Tc oc ns{T = 0)/m* 
is no longer valid there. The scaled curve Tc (x ns(T = 0)/m* crosses T* at 
optimal doping. Our result for Tc and T* already yield the observed behav- 
ior. We expect that a selfconsistent inclusion of phase fluctuations in FLEX 
will give better agreement with experiment, in particular a more rapid de- 
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Fig. 3. Doping dependence of the temperatures T* where Cooper pairs are 
formed, of where Cooper pairs becomes phase coherent, and of the phase 
stiffness ns{T = 0)/m*. T* and ns{T)/m* are calculated using FLEX and 
Tc is calculated using the Kosterlitz-Thouless criterion for the vanishing of 
the phase coherence. 

crease of Tc for x — > 0. The extended FLEX theory shown in figure I should 
give for the underdoped superconductors Tc oc Ug, ns{q=0,ih'n=0) = at 
Tc to obtain no Meissner effect above Tc, and the pseudogaps in the single- 
particle spectral density. Note, the coupling of the phase to its conjugate 
variable, the charge density, becomes important as the doping goes to zero. 
This may cause Tc —> for x —> 0. 
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